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Abstract 

We discuss the contribution coming from the infrared region to NLO matrix elements 
and/or coefficient functions of hard QCD processes. Strictly speaking, this contribution is 
not known theoretically, since it is beyond perturbative QCD. For DGLAP evolution all 
the infrared contributions are collected in the phenomenological input parton distribution 
functions (PDFs), at some relatively low scale Qo; functions which are obtained from a 
fit to the 'global' data. However dimensional regularization sometimes produces a non- 
zero result coming from the infrared region. Instead of this conventional regularization 
treatment, we argue that the proper procedure is to first subtract from the NLO matrix 
element the contribution already generated at the same order in a s by the LO DGLAP 
splitting function convoluted with the LO matrix element. This prescription eliminates 
the logarithmic infrared divergence, giving a well-defined result which is consistent with 
the original idea that everything below Qo is collected in the PDF input. We quantify 
the difference between the proposed treatment and the conventional approach using low- 
mass Drell-Yan production and deep inelastic electron-proton scattering as examples; and 
discuss the potential impact on the 'global' PDF analyses. We present arguments to show 
that the difference cannot be regarded as simply the use of an alternative factorization 
scheme. 



1 Introduction 



While studying the perturbative QCD (pQCD) description^ of the Drell-Yan production of low- 
mass pairs, in preparation for the interpretation of the forthcoming LHC measurement^ 
we unearthed a puzzle in the conventional procedure to remove the infrared divergencies. The 
puzzle is not just confined to the description of Drell-Yan production, but occurs for other QCD 
processes, such as deep inelastic electron-proton scattering, etc. The resolution of the puzzle 
can lead to numerical corrections at low scales. As such, it has the potential to influence the 
global PDF analyses [3]- jg]. 

To introduce the problem, recall that pQCD calculations are based on factorization theo- 
rems. The cross section of a hard process is calculated from the convolutions of parton distri- 
bution functions (PDFs) with the cross sections of the various hard subprocesses. The PDFs 
satisfy DGLAP evolution in InQ 2 starting from an initial input at some low starting scale Qq. 
The input values are known from 'global' fits to all the available deep inelastic and related hard 
scattering data. Below Qq the higher-order a s and higher-twist and other corrections, become 
too large for the use of pQCD. 

In the collinear approximation, at leading order (LO), each additional power of small a s is 
compensated by the large logarithm of a hard scale (such as InQ 2 ). At next-to-leading order 
(NLO) we have an a s term in the hard matrix element and/or the DGLAP splitting function 
unaccompanied by a large log. At NNLO the a 2 s terms do not have large logs, and so on. 

There is no problem with this procedure at LO — due to strong /^-ordering, when the 
DGLAP evolution is started from Qq, we do not encounter contributions from the kt < Qq 
domain. On the other hand, already inside the NLO cell there appears a contribution with 
kt < Qq- How to treat the (infrared) contribution from this large distance domain is the 
subject of this note. 

We discuss two alternative treatments of the infrared domain. First, we adopt a practical 
approach, which we shall call the 'physical' prescription, see Section [2} To avoid double count- 
ing, we have to subtract from the NLO diagram the a s *log contribution that has already been 
generated by LO DGLAP evolution. As we will show, this completely removes the infrared 
divergency from the matrix element. To be precise, after the subtraction, there remains a very 
small contribution coming from the k t < Qq domain of 0{(Q\j /J>%), where \ip is the factorization 
scale, which we will discuss in Section [2j 

The second treatment of the infrared divergency, which we will call the 'conventional' pre- 
scription, is based on dimensional regularisation, see Section [3] Here just the 1/e pole, which 
corresponds to the anomalous dimension generated by LO DGLAP evolution, is subtracted 
from the NLO contribution. In this case, after the subtraction, we are left with a contribution 
from the kt < Qq region which does not vanish as Qq — > 0. So, surprisingly, the two treatments 
of the infrared domain yield different results. 

X A preliminary study can be found in [1]. 

2 See, for example, the preliminary LHCb data [2J. 
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Sections [2] and [3] introduce the two treatments in more detail, using the NLO coefficient 
function of Drell-Yan production as an example. Then in Section 4 we quantify the difference 
of the contributions from the infrared domain taking again the NLO coefficient function in Drell- 
Yan production and, in addition, the 7*g coefficient function, C g , of deep inelastic scattering, 
as topical examples. Section [5] (as well as the latter part of Section 4) discusses the origin of 
the difference. 



2 Calculation at NLO — the physical approach 

Schematically, we may write the cross section in the form 

da/d 3 p = J dxidx 2 PDF(x 1 ,fi F ) \M(p; fM F , /Ar)| 2 PDF(x 2 , fi F ) , (1) 

where a sum over the various pairs of PDFs is implied. The matrix elements squared, |.M| 2 , 
describe the cross sections of the elementary partonic subprocesses. However, the problem is 
that using pi) we do not know the factorization scale /ip at which the PDFs are measured. 
Moreover, in the low x region, the PDFs strongly depends on the choice of p, F . 

In general, after the summation of all perturbative orders, the final result should not depend 
on the choice of [i F that is used to separate the incoming PDFs from the hard matrix element. 
Thus we need to account for the NLO, NNLO,.. .corrections. Contributions with low virtuality, 
q 2 < fi F , of the incoming partons are included in the PDFs, while those with q 2 > /j, F are 
assigned to the matrix element. 

Let us start with the LO expression for the cross section. In the collinear approach, the 
cross section has the form 

<r(ji F ) = PDF(/i F )®C LO ®PDF(/i F ). (2) 

The effect of varying the scale from m to fi F , in both the left and right PDFs, can be expressed, 
to first order in a s , as 

a{ji F ) = PDF(m) ® (c h ° + gin f^j (P lc{ tC LO + C LO P right ^j ® PDF(m), (3) 

where the splitting functions P r i g ht and Pi e & act on the right and left PDFs respectively. Recall 
that in calculating the a s correction in rt3J), the integral over the transverse momentum (vir- 
tuality) of the parton in the LO DGLAP evolution was approximated by the pure logarithmic 
dk 2 /k 2 form. That is, in the collinear approach, the Leading Log Approximation (LLA) is used. 

Now we turn to the expression for the cross section at NLO. First, we note that the original 
Feynman integrals corresponding to the NLO matrix element, C NLO , do not depend on fi F . 
However, we will see below how a scale dependence enters. At NLO we may write 

a(ji F ) = PDF(fi F )® (C L0 + a s C c N o ^ r )®PDF(/i F ), (4) 
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Figure 1: Diagrams (a,b) show NLO subprocesses for Drell-Yan production resulting from the 
splitting of the upper PDF (or 'right' PDF, in the notation of p])). 



where we include the NLO correction to the coefficient function. However to avoid double 
counting we must subtract the DGLAP-generated a s term of ([3]). This cancels the infrared sin- 
gularity, and gives a non-singular result which can be safely integrated in normal 4-dimensional 
space. 

2.1 Example: Drell-Yan production at NLO 

We illustrate the general procedure in terms of the specific example of Drell-Yan production 
of a /i + yU~ pair of mass M. At LO the subprocess is qq — > 7*, while at NLO we have the 
2—7-2 subprocesses qq — > (77* and gq — > 57* shown in diagrams (a,b) of Fig. [Tj These 
NLO contributions are now calculated with better, than LLA, accuracy. However part of these 
contributions are already included, to LLA accuracy, in the second term in ([3]), where the 
splitting functions 

-Pright = Pqq H~ Pqg j -Pleft = Pqq + Pqg (5) 

act on the right and left PDFs respectively So this LLA part should be subtracted from 

£<NLO 

For Drell-Yan production at low x, the majority of quarks/antiquarks are produced via the 
low-x gluon-to-quark splitting, P gg . That is, the most important NLO subprocess is gq — > q^*- 
The corresponding cross section, for the production of a pair of mass M and rapidity Y, 

has the schematic form 

= / dM* PDF^O PDF,(, 2 ) Si^s - M 2 ) S gh J - Y) , (6) 

where xj = Xi in the LO case, but xj are the true momentum fractions (x + and x~) carried by 
the photon at NLO. The subprocess cross sections have the form 



X ^-LO , -NLO fry\ 



We may equally well have incoming q's in P r i g ht and incoming g's in P\ 



left ■ 
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where, for the main NLO subprocess, we have 



da qg 


1 a 2 a s 


s t 2M 2 u 


dt 


9 s 2 


t s st 



Since M 2 = zs and u = M 2 — s — t, this gives 



da(gq — > q^*) a 2 a s z 1 



d\t\ 



9M 2 \t\ 



((1 -zf + z 2 ) + z< 



M 4 



2z< 



M 2 



(8) 



(9) 



In order to calculate the inclusive cross section da/dM 2 , it seems that we have to integrate 
over t starting from t = 0. If this were necessary, then we would face an infrared divergency. 



2.2 Treatment of the infrared region in the physical approach 



When dealing with PDFs we avoid the problems of confinement and interactions at large dis- 
tances. We start with some phenomenological input at a relatively large Qo, that is Ql ^> Aq CD , 
and consider just the evolution of the PDFs with Q 2 increasing from Qq to the factorization 
scale fi 2 F . Everything below Q is absorbed in the phenomenological input PDF. Accepting this 
logic, it is natural to replace the lower limit t = by the same Qq, and never to consider the 
contribution from low virtualities, k 2 < Qq. 

The subtraction of the LO DGLAP contribution from the NLO da/dt completely eliminates 
the 1/t singularity as t — > by introducing a 0(|t| — ftp) function in the first term of Let 
us explain, more explicitly, how this happens. We denote the gq — > qj* 0(a s ) contribution by 
The result of the explicit calculation of the gq — > q^* cross section is shown in d9J). On the 
other hand, this contribution may be written as the term generated by LO DGLAP evolution 
and the remaining NLO part, that is 



d*w = dd£? + <%> ® ^ 9 °dt (io) 



The last term reads 



where here DGLAP evolution has accounted for all virtualities \t\ = k 2 < fi 2 F ; and where the 
contribution of k 2 < Qq is hidden in the phenomenological input PDF. After the subtraction 
of this LO DGLAP generated term, the remaining contribution of ^ is 



da^Jf a 2 a s z 1 



d\t\ 9M 2 \t 
which has no singularity as t — > 0. 



\{l-z) 2 + z 2 ] 9(1*1-/4) + z 2 ^ - 2z 2 ^- 2 



(12) 



Recall that we need to subtract the DGLAP generated contribution of ( 11 ) in order to avoid 



the double counting. This subtraction is unambiguously defined. It is done in terms of t (and 
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not k t ) since the original DGLAP evolution was written as an evolution in k 2 = —t, and not as 
an evolution in kt- 

A small uncertainty is still present coming from the treatment of the non-singular terms 
in (12). We return to the choice of the lower limit of integration. It is not evident whether 
we should integrate from t = or from \t\ = Qq. From the ladder diagram, Fig. [2^a), which 
is of exactly the same form as that in LO DGLAP evolution in which everything below Qq is 
absorbed in the input, we are led to integrate from \t\ = Qq. However, the NLO hard subprocess 
cross section also includes non-ladder diagrams, Fig. |2|b) with an s-channel quark, for which 
the \t\ < Qq domain does not correspond to low virtuality. Thus the true uncertainty due to 
the non-singular terms is 0(Qq/ which may be neglected for fi F ^> Qq. 
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Figure 2: (a) Ladder and (b) non-ladder diagrams. 



3 Calculation at NLO — conventional approach 

We now consider the conventional prescription for the evaluation of 

(13) 



d&{gq — > <Z7*) a 2 a s z 1 \ rr . s0 , t 2 t 



>)*+*)+ 



d\t\ 9M 2 \t 

Recall, that to calculate the inclusive cross section da/dM 2 , we appear to have to integrate over 
t starting from t — 0, and that we face an infrared divergencjj^J In the conventional approach 
the integral is regularized by introducing the 4 + 2e dimensional space [91 [10l |TT] . Then the 
contribution from very small t produces a 1/e pole, which is absorbed into the incoming PDF. 
However, in 4 + 2e dimensional space, the number of gluon states with transverse polarisation is 
2 + 2e, giving e in the numerator of the matrix element. Besides this, there is an e dependence 
coming from the decomposition of the phase space factor, like (1 — z)~ e . So finally we have an 
e/e term, which produces a non-zero result as e — > 0. 



4 In perturbation theory this divergency is avoided by the (small) current quark mass, m q . However, in 
practice, the light quark mass is neglected. 
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Simultaneously, we have to consider the same diagram generated by LO DGLAP evolution, 
which gives a 1/e pole. The 1/e poles cancel, but we are left with the e/e term, which produces 
a non-zero result as e — > 0. Unlike the physical approach, this non-zero term does not vanish 
as 0(Ql/ fi 2 F ), and does not vanish for Ql <C /sp. 

Thus, suprisingly, we find that the physical and conventional prescriptions yield different 
results. We quantify this difference in the next Section, and then in Section [5] we discuss its 
origin. 

Sometimes the exact treatment of the infrared singularity does not matter. In particular, 
the emission of a soft gluon has a logarithmic infrared divergency which, according to the 
Bloch-Nordsieck theorem, is exactly cancelled between the real emission contribution and the 
virtual loop correction^} However, in NLO Drell-Yan or in the 7^ coefficient function for the 
structure function F2 of deep inelastic scattering, for example, the difference in the treatment 
of the infrared singularity does matter. 



4 Quantitative estimates 

To illustrate the numerical differences of the two treatments of the infrared contribution, we 
study, as examples, Drell-Yan production and the structure function F2 of deep inelastic scat- 
tering. 



4.1 Drell-Yan 

The Drell-Yan cross section calculated using the Vrap code [12j and MSTW08 NLO PDFs [3] 
is shown by the continuous curve in Fig. [3] for the LHC energy y/s = 7 TeV and a mass M = 6 
GeV of the produced pair. As usual, the factorization scale is chosen to be fip = M. The 

long-dashed line is the remaining contribution of the NLO gq — > q^* subprocess calculated in 
the framework of dimensional regularisation. We see that it turns out to be negative, since, for 



/ip = M, the contribution formally generated by LO evolution, (11), is larger than the whole 
0(a s ) cross section. 

On the other hand, if we were to adopt the 'physical' approach and subtract the explicit 
DGLAP-generated term (as discussed in Section 3), then we obtain the result shown by the 
short-dashed curve, which is significantly different to that obtained using dimensional regu- 
larisation. Moreover, we see the uncertainty coming from the \t\ < Qq domain is practically 
invisible in this figure — essentially identical results are obtained by integrating the last two 
(non-singular) terms in [...] in Q with a lower limit \t\ = Ql = 1 GeV 2 (dotted curve) and 
with t = (short-dashed curve). 



^This is the justification for the + prescription used for the 1/(1 — 2;) terms in the DGLAP splitting functions. 
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Figure 3: The LO+NLO cross section for Drell-Yan production of a pair of mass M = 6 

GeV and rapidity Y at pp collider energy of 7 TeV, obtained using MSTW NLO PDFs [3J. The 
two dashed curves correspond to the two different treatments of the infrared singularity in the NLO 
contribution. The lower dashed curve is obtained using 'conventional' e regularisation, whereas the 
upper dashed curve is obtained using the 'physical', approach working in normal 4 dimensional space. 
For the later we show the result for \t\ > (short-dashed) and imposing a cut \t\ > 1 GeV 2 (dotted) 
— the minute difference is invisible on the plot. 



4.2 Deep inelastic scattering 

An analogous situation is shown in Fig. [4] and [5] for the NLO gluon contribution to F 2 in DIS. 
The calculation is very similar to that for the Drell-Yan case since these two subprocesses are 
closely related by crossing symmetry. As before, the explicit subtraction of the term generated 
by the LO DGLAP evolution leads to a result for the 7*g coefficient function, 

C g = T R ([{l-z) 2 + z 2 ]\n- z + 6z(l-z) - lj , (14) 

that is different to that obtained by the conventional prescription, 

C g = T R [\{l-z) 2 + z 2 ]\i^^ + 8z(l-*) - 1^. (15) 

Here the notation of [H] has been used. The additional ln(l— z) singularity reflects the threshold 
singularity; as z — > 1 the phase space available for the quark and antiquark is proportional to 
(1 — z). In e regularization this threshold singularity produces the factor 

(l-z)- e = l-eln(l-z) + 0(e 2 ), (16) 



S 



which, on being multiplied by the 1/e pole, gives a constant ln(l — z) term, that does not vanish 
as e — > 0. Unlike the conventional treatment, the physical prescription has, at most, only a 
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Figure 4: The continuous curve is the LO+NLO result for the structure function F 2 (x,Q 2 ) for 
electron-proton DIS obtained from MSTW08 NLO PDFs [3J. The long-dashed and short-dashed 
curves are the NLO contributions arising from photon-gluon fusion (7*^7 — > qq) for the 'conventional' 
and the 'physical' treatment of the infrared 1/t divergence, respectively. 



very small 0(Qq/ ' [i 2 F ) infrared contribution coming from the non-singular terms in the infrared 
region. 

To gain insight into the difference between the 'conventional' and 'physical' treatments of 
the infrared divergency, we discuss the t-dependence of the j*g — > qq matrix element. Adopting 
the 'physical' procedure we start with the explicit perturbative QCD result and subtract the 
LO DGLAP-generated contribution. The result, after the photon polarizations have been 
convoluted with g^, may be written in the form 

\M\ 2 = \M t \ 2 + \M U \ 2 , (17) 

whereS 

\M t \ 2 = [{I - z) 2 + z 2 ]^Q(\t\ ~ ~ 1- (18) 
Similarly for Ai u with u interchanged with t. The result is non-singular as t — > 0. Numerically 



Integrating ( 18 ) over t 1 and adding the term which originates from Fl (where there is no infrared divergency), 
leads to (TiT). 
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Figure 5: The fractional correction to F 2 arising from using the 'physical', rather the 'conventional', 
treatment of the infrared divergency in the 7*g coefficient function, shown as function of Q 2 for two 
values of x, and for two different sets of NLO PDFs, namely MSTW08 [3] and CT10 [4]. 



the non-singular contribution coming from the < \t\ < Qq domain is not important^ After 
integration over t, this contribution vanishes as 0(Ql/fip) as Qo — > 0. Indeed, the 'physical' 
prescription eliminates the logarithmic infrared divergence, giving a well-defined result which 
is consistent with the original idea that everything below Qo is collected in the PDF input. 

On the other hand, the 1/t singularity in |.Mt| 2 , in the 'conventional' treatment, is first 
regularized by working in 4 + 2e space. After the subtraction of the analogous DGLAP-initiated 
contribution, the 1/e term is cancelled. However, the resulting oo — oo subtraction leaves some 
non-zero remaining contribution, which can be viewed as effective non-local vertices, which act 
as a S(t) contribution. Recall that analogous vertices arising from an ultra-violet divergence, 
play the role of local counter terms. Dimensional regularization provides us with the possibility 
to introduce such counter terms in a way which does not spoil the renormalizability of the 
theory. However, the analogous effective vertices of infrared-origin are unacceptable. They 
are contrary to the confinement property of QCD. Due to the confinement property it is not 
possible to continue the 1/t behaviour below some low scale of C(Aq CD ). Recall that in the 
global PDF analyses this low t contribution is included phenomenologically in the input PDFs 
at the starting scale Qo, which is taken to be larger than Aqcd- 

7 For the analogous Drell-Yan process we saw in Fig. [3] that the difference between the results obtained by 
integrating over the regions |£| > and |i| > Qq = 1 GeV 2 were invisible on the plot. The same insensitivity 
of the results, when using the 'physical' treatment of infrared domain, applies also to the short-dashed curve in 
the F 2 plot of Fig [4] 
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5 Discussion 



Our approach has some similarities with the 'subtraction method' [T3l [H] where to regularize 
the infrared divergency one subtracts from the real cross section of n + 1 parton production, 
da^ +1 , a term da A , which has exactly the same divergency. Simultaneously the da A term is 
added to the virtual loop correction for the n parton cross section, da^. After this both the 
real and the virtual contribution have no infrared divergency, and so may be calculated in the 
normal 4 dimensional space (without the e decomposition). 

In comparison with this 'subtraction method', where the subtraction term is introduced 
artificially, here we subtract the contribution really generated by DGLAP evolution. On one 
hand, this guarantees the absence of double counting, while on the other hand, it eliminates 
the infrared divergency. Note, also, that in our particular examples, like the NLO correction 
to i<2 due to photon-gluon fusion, or the NLO correction to the Drell-Yan cross section due to 
gq — > 97* subprocess, there is no analogous LO contribution caused by gluons. Therefore there 
is no NLO virtual correction, like da^ , in the 'subtraction method'. 

Let us return to the comparison of the two treatments of the infrared domain discussed 
in this paper: that is the 'conventional' e regularisation subtraction and the 'physical' direct 
DGLAP subtraction methods. Why do the results differ? The reason is subtle. In the 'physical' 
approach we do not need to assume that perturbative QCD is valid for virtualities k 2 < Qq, 
and, in particular, in the confinement domain. We evolve starting from some phenomenological 
'input' distribution at /i 2 = Qq. By considering different choices of Qq we verify that the 
contribution coming from the k 2 < Qq domain is small, and vanishes as Q^ — > 0. On the 
other hand, the conventional dimensional e regularisation procedure continues the 1/t singular 
behaviour below Aq CD , giving finally a non- vanishing infrared contribution, which may be 
viewed as a 8(t) term during the t integration. However, the part coming from low virtualities, 
k 2 < Qq, has no reason to have perturbative form, and, moreover, does not vanish as Ql — > 0. 

We elaborate this in a little more detail. Dimensional regularization is clearly a valuable 
tool in the ultraviolet region. But, here, we are concerned about its use in the infrared domain. 
The hope was that for infrared stable quantities (like F 2 ) the result will not depend on the 
contribution from the large distance domain; the real and virtual loop contributions will cancel 
each other. Sometimes this is true. In particular, due to the Bloch-Nordsieck theorem, we 
can use the so-called 'plus prescription' which eliminates the 1/(1 — 2) singularity in splitting 
functions. 

However, there are examples where 'conventional' dimensional regularization gives a non- 
zero infrared contribution of the form of e/e, which is the correct mathematical result of a 
particular perturbative QCD diagram calculation, but which makes no physical sense for real 
QCD with confinement. Here we have presented two examples (the NLO 7*g coefficient function 
in Drell-Yan production and C g in DIS) and have shown that the subtraction of the LO DGLAP- 
generated contribution completely eliminates the infrared singularity; the cancellation occurs 
between terms in C NLO and in the convolution the C h ° ® P L °. Recall that this subtraction 
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is necessary to avoid double counting. After this, the integral can be calculated in normal 
4-dimensional space, and, surprisingly, as demonstrated above, the results differ from those 
found in the 'conventional' approach. 

Similarly, using the conventional e regularisation procedure, the NLO splitting functions, 
in particular P^ LO (z), of DGLAP evolution may also contain contributions of infrared origin 
which are inconsistent with confinement. On the other hand, in the 'physical' approach, the 
absence of the infrared singularity is due to the cancellation of the contribution to the NLO 
splitting function arising from the direct calculation of the NLO Feynman diagram and the 
contribution already generated by LO DGLAP evolution. That is, it is due to cancellations of 
the contributions occurring in p NLO and in P LO <g> P L °. A more detailed discussion is given in 
the Appendix. The derivation of these new infrared-corrected NLO splitting functions is more 
complicated than the treatment of the coefficient functions that we have used as examples, but 
it is still possibl^J 

In summary, the problem cases for the 'conventional' treatment of the infrared divergence are 
the NLO coefficient and splitting functions, in particular the C g and P gg functions, which leave 
'erroneous' non-singular contributions. The main effect will be on the behaviour of the gluon at 
low scales, a domain where the present global parton analyses show unexpected 'valence-like' 
x behaviour of the gluon; see, for example, [SI SI IS] - In order to carry about such a 'global' 
investigation it will be necessary to compute, not only the coefficient functions, but also the 
NLO splitting functions using the 'physical' treatment of the infrared domain. 

One might think, say, in the DIS example that we considered, that the infrared corrections to 
the coefficient and splitting functions could compensate each other, and that the transition from 
the 'conventional' coefficient function (which includes a non-vanishing infrared contribution 
inconsistent with colour confinement) to that obtained in a 'consistent' treatment, may be 
effectively accounted for by a re-definition of the PDFs, and thus considered as an alternative 
'factorization scheme'. However such an alternative scheme, which may correspond to new 
coefficient functions, prescribes a new set of splitting functions to describe the evolution of the 
re-defined PDFs, which are, in general, different to the infrared-corrected splitting functions. 

Even without explicit knowledge of the NLO splitting functions corresponding to the 'phys- 
ical' treatment of the infrared domain, there are several ways to see that compensation is not 
possible and that the results cannot be regarded as simply a scheme change. A hint, that this 
is so, is the presence of the additional ln(l — z) singularity in the 'conventional' C g of (15). 
An argument is to consider the non-singlet channel. Here, due to flavour conservation, the 
PDF normalization is fixed, and we have no possibility to re-define the non-singlet PDF via 
the admixture with some gluon or other quark flavour contribution. In other words, there is 
no chance to compensate the correction to NLO coefficient function by the re-definition of the 
(non-singlet) parton density and/or by the correction to the NLO splitting function, which is 
exponentiated, that is iterated many times already at NLO, especially at small x. It is not a 
NNLO effect. 



3 We are at present tackling this problem. 
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Alternatively, to be more explicit, we may consider the quark contribution to the gg splitting 
function. This is very similar to the coefficient function of 7* g fusion. The part of the splitting 
function of (inconsistent) infrared origin gives a contribution, A7, to the anomalous dimensions, 
7, - - it therefore exponentiates during the evolution, which is especially long in the low x 
region. In terms of the anomalous dimensions, the inconsistency A7 is strongly enhanced in 
the low x region by the small value of the analytic continuation of the Mellin moment u, recall 
A7 ~ const. /w. Compensation is clearly not possible over a long interval of Q 2 . 

Another good, related, example is evolution in QED (where there is no confinement) and 
where the electron density is well defined. Here the evolution was first considered by Gribov 
and Lipatov [J5] in terms of the Leading Log Approximation (rather than the 'factorization 
theorem' approach), which allows one to trace each step of the evolution and the contribution 
of each Feynman diagram explicitly. As there is no confinement, the electron and photon fields 
are well defined. Thus the only possibility to re-define the 'partons' is a field renormalization 
Z(q 2 ) factor which does not depend on the momentum fraction z. Of course, one may consider 
some mixture of electron and photon PDFs, but we discuss the original electron and photons. 
For these distributions the only way to avoid double counting is the 'physical' approach. We 
emphasize that it is useful to first consider the problem of double counting in QED where 
everything is well defined. After this it is straightforward to discuss what happens in QCD. 

A more general observation is that in our approach we work with partons defined by OPE 
operators, and the corresponding Feynman diagrams. For these quarks and gluons we evaluate 
the Feynman diagrams in both the 'conventional' and the 'physical' approach. In either case, 
we first have to calculate the splitting and the coefficient functions for the normal (OPE) quarks 
and gluons (for which we have Feynman rules and where each line in the diagrams has definite 
quantum numbers). Only after this may we define the new q' and g', and new splitting (and 
coefficient) functions corresponding to an alternative factorization scheme. So differences in 
the infrared treatments cannot be attributed to a scheme change. 

We conclude that differences in parton behaviour that follow from an analysis using a 
'physical' treatment of the infrared domain cannot be reproduced by a factorization scheme 
transition from the conventional PDF analyses. (The differences are expected to be mainly in 
the gluon at low scales.) The only way to account for the appropriate treatment of the infrared 
region is to perform a new global analysis using a complete set of corrected coefficient and 
splitting functions. 

Appendix: Infrared contributions to NLO splitting func- 
tions 

In the axial gauge the only infrared singularity of a splitting function is due to the ladder-type 
(box) diagram (plus a self-energy contribution oc 5(1 — z) dzP(z) which can be determined 
from this diagram in usual way based on the flavour and energy- momentum conservation laws). 
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The singularity has a logarithmic form dk 2 jk 2 and is exactly equal to the contribution P LO cg)P LO 
generated by LO DGLAP evolution, since there is no other dimensionful parameter apart from 
Q 2 . To avoid double counting, we have to subtract the LO DGLAP-generated contribution 
from the NLO splitting p NL °. This subtraction eliminates the infrared singularity of P NL °. 
However the remaining part may still contain some non-singular contribution of infrared (large 
distance) origin. 

Recall that DGLAP evolution is performed starting from some relatively large virtuality 
Ql > Aq CD . Everything below Q is collected in the phenomenological 'input' parton distribu- 
tions, PDF(Qo)- It would be best to calculate the NLO box diagram also starting the integral 
over k 2 from Q^; that is to consider dk 2 /k 2 . On the other hand, in such a case we will in- 
troduce a new dimensionful parameter which will strongly complicate the logarithmic DGLAP 
evolution. If we assume that the starting scale is much smaller than the final factorization scale 
li 2 F = Q 2 , that is Qq <C Q 2 , then we may consider the non- singular contribution from the region 
with k 2 < Qq as a power, 0(Ql/Q 2 ), correction, and neglect it, together with the other power 
corrections. That is to consider the J Q dk 2 /k 2 . 

The problem is that within the conventional approach we deal with singular integrands. 
Using the dimensional, e, regularization, this results in a 1/e pole plus some constant terms 
(with respect to e) of infrared, e/e, origin. Finally we subtract the 'anomalous dimension' 
function]^] generated by the LO DGLAP evolution; that is we eliminate the 1/e pole. However, 
this prescription does not cancel the constant terms (with respect to e) of infrared, e/e, origin. 

If, on the other hand, we adopt the 'physical' approach and integrate the non-singular 
expression given by the box diagram minus the DGLAP-generated P LO <S> P LO contribution 
(or put some infrared cutoff go on the integral J^dtf/k 2 ) we will never obtain such terms 
of infrared origin, even for go — > 0. We thus conclude that the above LO DGLAP-generated 
subtraction will give the correct NLO splitting functions, and that the term coming from very 
large distances (greater than C(1/Aqcd)) in the 'conventional' procedure should not be there. 
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